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Abstract 

We study the scenario of electroweak baryogenesis mediated by nonsuperconducting 
cosmic strings. This idea relies upon electroweak symmetry being restored in a region 
around the core of the topological defect so that, within this region, the rate of baryon 
number violation is enhanced. We compute numerically how effectively baryon number 
is violated along a cosmic string, at an epoch when the baryon number violation rate 
elsewhere is negligible. We show that _B-violation along nonsuperconducting strings 
is quite inefficient. When proper accounting is taken of the velocity dependence of 
the baryon number production by strings, it proves too small to explain the observed 
abundance by at least ten orders of magnitude, whether the strings are in the friction 
dominated or the scaling regime. 
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1 Introduction 



Considerations about the formation of light element abundances when the temperature of 
the Universe was about 1 MeV lead us to conclude that the difference between the number 
density of baryons and that of antibaryons is about 10 -10 if normalized to the entropy density 
of the Universe 

The theories that explain how to produce such a tiny number go generically under the 
name of theories of baryogenesis, and they represent, perhaps, the best example of the 
interplay between particle physics and cosmology. Until now, many mechanisms for the 
generation of the baryon asymmetry have been proposed. Grand Unified Theories unify the 
strong and the electroweak interactions and predict baryon number violation at the tree 
level. They are, therefore, good candidates for a theory of baryogenesis. There, the out- 
of-equilibrium decay of superheavy particles can explain the observed baryon asymmetry 

§• 

An alternative possibility is represented by electroweak baryogenesis, which has been 
the subject of intense activity in the last few years || ||. The Standard Model (SM) itself 
provides one of the key ingredients for baryogenesis || : violation of baryon number via the 
anomaly |J , which also constrains some models in which the baryon asymmetry is generated 
at a very high energy scale. 

The underlying idea is this. If the electroweak phase transition is of the first order, 
the minimum of the Higgs potential associated with the symmetric phase is separated from 
the local minimum of the broken phase by an energy barrier. At the critical temperature T c 
both phases are degenerate in energy and at later times the broken phase becomes the global 
minimum of the potential. The universe supercools in the symmetric phase, and the phase 
transition then proceeds by nucleation of critical bubbles. As the bubble walls separating 
the broken from the unbroken phase pass each point in space, the order parameter changes 
rapidly, leading to a significant departure from thermal equilibrium. CP-violating sources 
may be locally induced by the passage of the bubble wall and baryogenesis is enhanced if 
the CP-violating charges can efficiently diffuse in front of the advancing bubble wall, where 
anomalous electroweak baryon number violating processes induced by sphaleron transitions 
are unsuppressed 0. 

In the SM, the electroweak phase transition is too weakly first order to assure the preser- 
vation of the generated baryon asymmetry, as perturbative and non-perturbative analyses 
have shown ||. For this reason, electroweak baryogenesis requires new physics at the elec- 
troweak scale. Low-energy supersymmetry is a well-motivated possibility, and there is still 
some room for a suitably strong phase transition in this scenario |J and for the production 
of the baryon asymmetry H. 

In this paper we will focus on a different mechanism to implement electroweak baryoge- 
nesis, in which the third Sakharov condition, loss of thermal equilibrium 0, is fulfilled by 
the evolution of a network of topological defects ||, H| . The advantage of this idea is that it 



does not depend upon the order of the electroweak phase transition, although the presence 
of stable topological defects still requires new physics beyond the SM. 

Topological defects are regions of trapped energy density which may be left over after a 
cosmological phase transition if the topology of the vacuum of the theory is nontrivial. In 
particular, cosmic strings are solitonic solutions in spontaneously broken field theories for 
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which there exist non-contractible loops in the vacuum manifold |TT . 

An essential ingredient of cosmic string mediated baryogenesis is the restoration of elec- 
troweak symmetry in a region surrounding the core of the string. This obviously requires 
some coupling of the fields responsible for the strings, which are presumed to have formed 
at a scale above the weak scale v, to the Higgs field. For example, the strings may arise 
from the breaking of an extra gauged U(l)' symmetry by the formation of a condensate of 
some other scalar field S; if the SM Higgs is charged under U(l)', electroweak symmetry is 
restored in a region about the string because this reduces Higgs field gradient energies within 



the string [12 



After the electroweak phase transition, baryon number violation may be efficient only in 
the symmetry restored region around the string, and not outside. If, in addition, there are 
CP violating interactions between the string and the background plasma, the motion of the 
string can disturb equilibrium in a way which creates a localized CP asymmetry. The CP 
asymmetry biases electroweak sphalerons to generate a baryon asymmetry [[|, in a manner 
analogous to the physics at a bubble wall in the conventional scenario with a first order phase 
transition. We anticipate that the asymmetry generated by cosmic strings will be smaller 
because only a fraction of space is swept out by the moving network of strings, and because 
the opposite walls of the string produce asymmetries of opposite sign that tend to cancel 
unless baryon number violation is very rapid inside the string. 

From these considerations one sees that a crucial requirement for string mediated elec- 
troweak baryogenesis is that the baryon number violating sphaleron transitions are fast in 
the regions near the strings. Further, this must be true for some part of the epoch where 
baryon number violation is very inefficient in the bulk plasma, because when baryon number 
violation is at all efficient in the bulk, any asymmetry produced by strings is subsequently 
erased. Hence, the mechanism relies on there being some temperature range in which both 

1. Sphaleron configurations large enough to be unsuppressed fit inside the region of elec- 
troweak symmetry restoration, so the baryon number violation rate is unsuppressed 
close to string cores, and 

2. sphalerons are suppressed in the bulk plasma, so baryon number produced by strings 
is preserved. 

To see what is needed to satisfy the above requirements, let us consider the typical length 
scale R s of symmetry restoration around the center of nonsuperconducting strings: 

1 



X 1 / 2 v(T) ' v ' 

where A is the quartic Higgs coupling and v(T) is the temperature dependent vacuum expec- 
tation value of the Higgs field. On the other hand, the sphaleron size in the unbroken phase 
is R sp h ~ (g 2 ^) -1 - Conditions 1. and 2. are then R sp h ^ R s and v(T)/T > 1, respectively. 
Together these imply 

X<g\ (2) 

which indicates that one needs a Higgs boson which is far too light and already ruled out by 
present LEP data which require A ~ g 2 [|T3[1 . It might be possible to satisfy this condition in 
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extensions of the SM; however such a weak scalar self-coupling leads to a very strong first 
order electroweak phase transition, in which case baryogenesis could proceed without strings 
anyway. 

It is therefore clear that a more rigorous treatment is needed to understand whether string 
mediated electroweak baryogenesis can really explain the observed baryon asymmetry. To 
answer this question, we have computed either the sphaleron energy or the sphaleron rate 
in the background of each sort of nonsuperconducting strings. For sphaleron energies we 



use techniques similar to those in [|J] and for sphaleron rates we use the techniques of [fL5 



The case of superconducting strings is more complicated, because it involves determining the 
size of currents on such strings in a realistic cosmological setting. Since we are not aware of 
any reliable estimates of the current carried by a superconducting string in a realistic early 
universe cosmological setting we have not attempted this here and we leave it for future 
investigation. 

Our findings indicate that the baryon number violation along a string is generally too 
slow for significant baryogenesis. This is because the region where electroweak symmetry 
is restored is never wide enough to allow symmetric phase sphaleron-like events. At best it 
brings the sphaleron rate to 1/30 of what it would be if the core region were of the order of 
the natural nonperturbative length scale 1/ctwT; and this only occurs in a rather special case 
where the Higgs field is charged under an extra U(l)', with a charge incommensurate with 
the charge of the scalar field which breaks the extra U(l)'. If the charge is commensurate, 
or if the Higgs field only interacts with the string generating fields via potential terms, then 
the sphaleron rate is insignificant. 

Moreover, if the network is in the scaling regime, the baryon number violation occurs in 
too small a total volume to explain the observed baryon asymmetry. If the network evolution 
is friction dominated, the network will be denser, but strings will move more slowly. As we 
discuss in some detail here, the baryon number production is suppressed by the second power 
of the string velocity, so this case also leads to insignificant baryogenesis. In either case we 
find that the generated baryon number is at least ten orders of magnitude smaller than the 
observed abundance. 

The paper is organized as follows: section 2 is dedicated to the computation of the baryon 
number violation on the strings for different cases; section 3 deals with the generation of CP 
asymmetry inside the strings; and section 4 with the efficiency of baryogenesis by a string 
network. Finally, in section 5 we present our conclusions. 



2 Baryon number violation on cosmic strings 

To determine how efficiently cosmic strings can generate the baryon number asymmetry, we 
need to know the rate of B violation by sphalerons inside the string, at times when the 
Higgs field condensate is large enough so B violation is negligibly slow outside the string. 
We characterize the efficiency of baryon number violation on the string by the Chern-Simons 
number (Nqs) diffusion constant per unit string length L, 

s , im <w*m - N Csi om 

t^oo Lt y ' 
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This is related to baryon number production by the string through 



L dt ~ 2 i T ' U 

where Np = 3 is the number of families and \x\ is half the difference in chemical potential 
between particle and antiparticle, averaged over the region of symmetry restoration. The 
sum runs over each doublet which couples to SU{2) L . We will often write Ti in terms of a 
dimensionless quantity Ki, 

T t = K ia 2 w T 2 . (5) 

For a very thick string, k\ would numerically be approximately o? w T 2 times the cross sectional 
area, although parametrically there is actually an extra power of a w [JlSj, |17], [18], [19], [2(J . We 
would say that baryon number violation is unsuppressed on the string if acj ~ 1. 

Since the only baryons which survive are those created after the Higgs VEV exceeds the 
temperature, we will focus on temperatures such that 4>q/T > 1, where (f> is the amplitude 
of the Higgs field condensate, (fil = 2$t$ (broken phase). We expect (and can check) that at 
lower temperatures, baryon number violation becomes less efficient along a cosmic string. For 
instance, a cosmic string today, at temperatures < lMeV, will not catalyze baryon number 
destruction except by extremely inefficient instanton processes, unless the region of symmetry 
restoration is very large-either > 1/awT, so sphaleron processes will be unsuppressed inside, 
or > 1/Aew; the scale at which the electroweak coupling becomes large and instantons 
become efficient. (This scale is enormous, 1/Aew > 1 meter.) 

Furthermore, we will focus on the case m# ~ mw, so the scalar self-coupling is of order 
the gauge coupling. Smaller values of m# are experimentally excluded, and larger ones can 
only decrease the efficiency of baryogenesis. In parametric estimates we will always write g, 
meaning either g, ^ g 2 + g' 2 , or vA. 

Except perhaps for superconducting cosmic strings, the width of the region of symmetry 
restoration on a cosmic string is of order the Higgs field correlation length l/g0o, which 
by no coincidence is the characteristic size of a sphaleron in the broken phase. So it is 
not a priori clear that sphaleron processes are unsuppressed on cosmic strings. Rather, we 
should check whether this is the case. We now endeavor to do so, case by case for various 
symmetry-restoring types of strings. 



2.1 Strings which force <£> = on the string core 

We begin with strings which force the Higgs condensate to zero along the core of the string, 
assumed much thinner than l/g<fto, but do not influence electroweak physics outside the 
string core. There are two obvious examples of such strings. 

In the first, the standard model is extended by a complex scalar field S, which is either 
a singlet under all gauge interactions, or transforms under an extra U(l)' symmetry for 
which all the standard model particles are uncharged. The effective potential for S and the 
standard model Higgs field $ is 

V{3, $) = A s (S*S - S 2 ) 2 - i(S*S - S 2 ){&<5> - $o) + A ($ f $ - $ 2 ) 2 , (6) 
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where we have assumed that the coefficient of the interaction term between the S and Higgs 
fields is negative; if it is not, electroweak symmetry will not be restored. Here and throughout 
we will use the notation in which capital letters S and $ denote complex fields, while lower 
case s and <fi are real components in the direction of the condensate, normalized such that 



*8 = 2^. 



Strings will form at some high temperature T ~ \J\ s Sq where the S field condenses. 
If the S field has no gauge interactions, the strings are global; if it transforms under a 
local extra U(l)' then they are local (Abelian Higgs model) strings. In either case, by the 



electroweak scale, the S field carries a condensate So = y\(S)\ 2 . The doublet Higgs mass 
parameter is — 2A$g + j(Sq — \S\ 2 ). In the bulk, where S = S , electroweak symmetry will 
break because the Higgs field mass term has the usual form, — 2A$q. However, in the core of 
a cosmic string where the S field condensate vanishes, the Higgs field feels a large positive 
mass squared ~ 7S5 2 , which may be enough to force $ = inside the string core. Whether 
it does so or not is determined by the competition between the Higgs field potential and 
gradient energies. The outcome of this competition depends on the values of the quartic 
couplings and So, which also determine the thickness of the string core. We will discuss this 
further, in particular the possibility that the S field string is fat, in the subsection on fat 
strings below. In the present section we will restrict our attention to the case where the S 
field string is thin compared to the inverse weak scale, 1/0q, and the string does successfully 
force $ = in its core. 



Another type of string which forces $ = in its core was considered in ref. ||12|| . The 
idea is that there is a new U(l)' gauge symmetry, broken by a scalar S with charge q s under 
£7(1)' but neutral under the SM gauge group. The SM Higgs field is assumed to also have 
a charge q<$, ^ under U(l)'. The U{1)' breaks at some scale above the electroweak scale, 
where the S field condenses with So 3> $o- A cosmic string will carry a U(l)' magnetic flux 
of 2n/q s , and the phase of the Higgs field gets shifted on parallel transportation around the 
string by 27rq<p/q s . For now we will assume that q^jq s is an integer. If it is not an integer, 
the Higgs field is forced to be small in a wider region around the string, and a Z magnetic 
flux will also be established; we treat this special situation in the next section. 

In the present case, the Higgs field minimizes gradient energies by having its phase wind 
by —l-Kq^jqs in going around the string, to compensate for the phase induced by the con- 
nection. This prevents azimuthal gradient energies in the Higgs condensate from occurring 
outside the string core. But inside, where a loop does not enclose the complete U(l)' magnetic 
flux, a nonzero Higgs condensate, if it existed, would possess a large, energetically expensive 
azimuthal covariant gradient, so the Higgs condensate is forced to zero. Writing the U(l)' 
field strength as Be = b(r)/(q s r) and the Higgs field as $(r, 6) = /i(r)$ exp(— iOq^/qg), the 
Higgs field gradient energy per unit length of string is 

= 2, f rdr (*§(!- WMjl + mr) ,\ (7) 
dL J [ qj r 2 J 

The very large coefficient on the first term at small r, where b ^ 1, forces h ~ here, but h is 
free to rise away from the core. Far outside the cosmic string, we can redefine the Higgs field 
to include the phase which compensates for the U(l)' connection, and electroweak physics 
then looks normal, even though $ is pinned to zero inside the core of the string. 
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The above scenario is phenomenologically constrained. Since the Higgs field is charged 
under the extra U(l)', and it couples to the standard model fermions, they must carry 
U(l)' charges as well . Because the Higgs field Yukawa couplings to fermions mix the 
generations, these charges must be generation- independent. Then the only anomaly- free 
charge assignment for the fermions is times their hypercharges. Thus the U(l)' gauge 
field is a standard-model-like Z', which is experimentally excluded for masses up to about 
700 GeV. The symmetry breaking scales of U(l)' and the electroweak sector have to be 
separated by at least one order of magnitude, which ensures that this string generation 
mechanism will always fall into the "thin string" case. (In the case of a U(l)' whose generator 
is not proportional to that of hypercharge, anomalies might be cancelled by adding additional 
heavy fermions. The mass limits on these other kinds of Z' bosons are similar.) 

For either of the two possibilities described above, the SM Higgs field is forced to be 
exactly zero inside the core of a narrow string of thickness r , and it rises to its normal 
condensate value (at the given temperature) outside. The Higgs field profile is that which 
minimizes the energy 

— = 2vr / rdr {& ti 2 (r) + V(h%)\ . (8) 

If we assume that dh/dr < h/r then the first term behaves as dr/r and the second 
behaves as rdr. Hence the potential term dominates at large r but at small r it is negligible. 
To see how the profile depends on r (the size of region where 5 C So, which is the width 
of the string), assume that the role of the potential is to force the Higgs field to reach its 
condensate value within a radius of r\ ~ l/g<&o, and can be neglected for r < r%. Then the 
energy is roughly 

dE rmn / dh \ 2 

*o / 1ruTZ\ rf ( lnr )- (9) 



2ndL 7inr \G?(lnr 

Each equal logarithmic interval of r contributes a comparable amount to the rise of <3>. Hence 
the size of [h(r) — 1] near r = n will depend on ro through 

K r )- 1(X n K 7T7TT- ( 10 ) 

m(n -r Q ) + 0(1) 

The smaller the string core is, the more of the Higgs field's rise occurs at small r, and the 
weaker the symmetry restoration at r ~ 1/g^Q. But the dependence is only logarithmic in 

At distances of order l/g<&o the suppression of $ below $o is incomplete; rather than 
$ ~ 0, we find $/$o ~ 1 — 1/ log( r / r o)- But r ~ l/g&o is the physical scale important for 
setting the sphaleron rate. This suggests that for ro much less than l/g&o, the energy of a 
sphaleron on a string will not be significantly less than in the absence of a string. 

To test this, we should determine the energy of a sphaleron in the string background. A 
sphaleron in the broken phase has a high degree of symmetry; a sphaleron Ansatz with gauge 
and Higgs field profile functions depending on radius alone ]2T| describes the saddle point 
unless A/ g 2 ^> 1 . However the string breaks the spherical symmetry down to cylindrical 
symmetry. A solution based on an Ansatz would need profile functions dependent on r and 
z separately. We prefer to find the sphaleron solution on the lattice. One writes down a 3D 
lattice Hamiltonian and seeks the saddle point solution corresponding to the sphaleron. We 
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a x g(f) 




vC9 


V0.4 


B 


1.24 


1.30 


1.37 



Table 1: Dependence of sphaleron energy E sph 
core string. 



47rB<j) Q /g, on lattice spacing for the thin 



use an 0(a 2 )-improved lattice Hamiltonian (a is the lattice spacing) so lattice spacing errors 
will begin at 0(a 4 ). 

We make a starting guess for a sphaleron configuration by the same technique used in 
|T5|| . Dissipative cooling will bring this configuration close to the sphaleron, but since the 



sphaleron is a saddle point and there are always round-off errors and lattice artifacts, the 
configuration will miss the sphaleron and cool to the vacuum. An algorithm to cool toward 



saddle points was developed especially for dealing with sphalerons on the lattice in |L4j] , but 
its use becomes cumbersome for improved actions. We use a new algorithm for cooling to 
saddle points, presented in Appendix |A]. We confirm that the candidate numerical solution 
for the sphaleron really is one, by measuring its Chern-Simons number using the technique 
developed in |23j. 

First, we confirm that the algorithm is working correctly by determining the sphaleron 
energy in the bulk, for X/g 2 = 1/8. On parametric grounds the sphaleron energy can be 
written as E sph = AnBcpQ/g, with B a pure number. We find B = 1.815 with good lattice 
spacing independence, in excellent agreement with the result using the sphaleron Ansatz, 
B = 1.82. 

We include the string by pinning the Higgs field to zero along a line of lattice sites. The 
value of B we obtain at three values of the lattice spacing are given in Table |I[ By "sphaleron 
energy" we mean the energy of a sphaleron in the string background, minus the energy of the 
string background alone. The weak lattice spacing dependence occurs because the effective 
core size is getting narrower, since our procedure makes it one lattice spacing wide. It is 
worth commenting that the real space distribution of the magnetic energy associated with 
the sphaleron remains fairly close to spherically symmetric; the sphaleron does not stretch 
out along the string significantly. This is also true if we force $ = by hand in some wider 
cylindrical region. 

These results suggest that the best hope for electroweak symmetry restoration is to have 
a fat string, which in turn requires symmetry breaking of the singlet field very close to the 
electroweak scale. This case demands a more careful analysis because the back reaction of 
electroweak fields on the fields generating the string cannot be neglected, and we will return 
to it in a later section. As for the case where the Higgs field is charged under an extra U(l)', 
we note that none of our lattices had as narrow a core as is actually required by the Z' mass 
bound, so the sphaleron energy is even higher than what we have found. 

The conclusion is that E sphy with string > 0.7£ sph , broken phase, getting closer to the broken 
phase value as the string's core becomes narrower. The relation between the sphaleron 
energy and the baryon number violation rate involves zero-mode contributions which will 
also be smaller for the sphaleron on the string, since it has less symmetry. In particular, 
the translational zero modes previously gave / d 3 x, bounded by the region of consideration. 
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This will become l 2 L, with L the length of the string and I given roughly by the distance the 
sphaleron can be translated away from being centered on the string, at an energy cost of T. 
We expect, roughly, l~ 2 ~ (gfo) 2 (£ sphi broken phase - -E sph , with string) /T. This estimate follows 
from assuming that displacing the center of the sphaleron slightly away from the string core 
raises the sphaleron energy by an amount quadratic in the displacement, and a displacement 
by the width of the sphaleron, ~ l/g<f) , changes the energy to the broken phase (no string) 
value. 

If we approximate the rotational zero mode contributions as being the same as for a 



sphaleron in the broken phase, then we may take the values determined in [24]. This gives 



K\ ~ exp(— 12) for v/T = 1. In fact the rate is smaller because the rotational symmetry is 
also broken, so the rotational zero modes will be partially removed. We made a preliminary 
calculation using the nonperturbative techniques of [[^|; the result was Ki ~ 10~ 6±0 ' 5 , which 
is uninterestingly small, so we did not refine the calculation further. 

We conclude that if a string only influences electroweak fields in its core, and the core 
width is less than the weak scale m^f, then it fails completely to allow efficient baryon 
number violation. 



2.2 Strings of non-integer magnetic flux 

In the previous subsection we saw that if the Higgs field transforms under an extra U(l)' 
with charge q^, and a string carries a U(l)' magnetic flux of 2ir/q s with q^jq s an integer, 
then the Higgs field is pinned to zero in the core of the string. The case where q$jq s is not an 
integer is much more interesting, and has been explored by Davis and Perkins ||12|| . Outside 
of the string, if there are no other gauge field condensates, the Higgs field phase changes by 
2nq ( f ) /q s on parallel transport around the string, using the gauge field as the connection. If 
the Higgs field condensate has a phase which changes by —2nn around the string, so that 
$(r, 9) = $ h(r) exp(—ni9) , and if n ^ q^/qs, as will be the case if q^/q a is not an integer, 
then there will still be a gradient in the azimuthal direction, leading to a gradient energy of 



2^1 J rdr L' 2 + (n - q*/q s ) 2 ^ } • (11) 



Since h(r) must eventually reach its asymptotic value of 1 to avoid an extensive potential 
energy cost, this leads to a logarithmically divergent gradient energy, the logarithm arising 
from J(l/r 2 )rdr. This can only be removed by having a nonzero Z boson field, with a net Z 
field magnetic flux of 2ir(n — q$,/q s )(2/g), to compensate. Since the total Z magnetic flux is 
fixed by the requirement that the Higgs gradient energy vanishes at large distances, / B z d 2 x 
is fixed. But the energy in the Z magnetic field is J B"^d 2 x, which is minimized (for fixed 
/ Bzd 2 x) by spreading Bz out over a large region of space. This leads to a larger region 
of symmetry restoration than in the previous section, and presumably a smaller sphaleron 
energy. 

The string is widest when q^jq s = 0.5 modulo 1. For this case, writing the Z field as 

* = (12) 
gr 
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higgs in string 
higgs in sphaleron _ 
gauge in sphaleron _ 
Z in string 



2 4 6 8 

radiusxg0 Q 

Figure 1: Higgs and gauge field profile functions for the string and the sphaleron. The 
Higgs field rises faster in the string, but the Z field rises more slowly than the gauge field of 
the sphaleron. In each case, A = g 2 /S. 

(approximating that the weak hypercharge gauge coupling g' is zero henceforth), the energy 
per unit length of string is 

dE r°° , f 1 z' 2 2 ,, 2 (z-l) 2 $ 2 h 2 x 4/l2 x2 

The solution is independent of the core radius of the string provided that it is small enough. 
We have therefore taken the r — > limit in this equation. 

We can solve for the profile functions z(r) and h(r) which minimize this energy, and 
compare them to the gauge and Higgs profile functions of a sphaleron; see Figure [I]. The 
figure is slightly misleading for two reasons: first, the string profiles would rise somewhat 
faster had we not made the approximation g' — 0; and second, the string profiles are for 
the radial direction in cylindrical coordinates, while the sphaleron profiles are for the radial 
direction in spherical coordinates. Hence more of the sphaleron lies within the string than 
the profile functions suggest. The string and the sphaleron are about the same size, and we 
should expect an order unity shift in the sphaleron energy in the presence of the string. The 
situation here is much more favorable than in the last subsection. 

However, the above argument is too naive in the present situation. Unlike the case in 
the last subsection, where the effect of the string could be treated as just a modification of 
the Higgs potential along the string core (shifting the mass up enough to pin $ = 0), here 
the string background already contains large gauge magnetic fields and azimuthal Higgs field 
gradients, which will influence the fields of the sphaleron configuration, potentially changing 




10 



the sphaleron energy substantially. Directly measuring the sphaleron energy shows that in 
fact it is significantly lowered. 

When the sphaleron energy becomes small (of order the temperature), the perturbative 
computation of Ti by the saddle point method becomes less reliable, and we are better off 
using a fully nonperturbative technique. If the rate is not too small, we can watch sphalerons 
occur during a real time evolution by tracking Chern-Simons number topologically ||15|| . If the 
rate is so small that no sphalerons will occur during a reasonable amount of Hamiltonian time, 
we can use a new technique based on finding the sphaleron free energy nonperturbatively on 



the lattice [p3| . It turns out in our case that for parameters which give a Higgs condensate 
of <fio/T = 1, the sphaleron rate is large enough on the string that the real time approach 
can be applied. 

We work on a lattice of spacing a = 2/(5g 2 T), using (O(a)-improved p5| ) scalar self- 
coupling A = 0.125g 2 and a value of the Higgs mass squared parameter which gives v/T = 1 
(technically, 2 = 2.5g 2 T 2 ). We set the Weinberg angle Q\y = 0, that is, we ignore weak 
hypercharge, so that we only have to deal with one U(l) gauge field. The volume is a 
3D L x x Lx x L 2 torus. We use L x = L 2 = 12.8/ g 2 T (a 32 3 lattice) and L 2 = UA/g 2 T, 
L x = 16/g 2 T (a 36 x 40 2 lattice). Both volumes are much larger than it takes to get 
continuum-like behavior in sphaleron rates P6[ , but we need to check whether the string's 
presence changes this picture. 

We fix a U(l)' background field appropriate for a half unit flux along a vertical line at 
(L/2, L/2, z). For topological reasons the integral of magnetic flux across the cross section 
of the box must equal zero, so there must be a return flux somewhere; we put it all on the 
plane x = (which is identified with x = L) or the plane y = (identified with y = L) in 
such a way that the connection is undisturbed from the string solution except along these 
planes. Using two box sizes allows us to check that the return flux is not affecting the results. 
Besides this fixed U(l) background, the approach is the same as used in fll5| . In particular 
we do not introduce extra degrees of freedom to generate hard thermal loop effects beyond 
those induced by UV lattice modes, although the technology exists [[0|]. This makes our 
results an overestimate, by of order a factor of 2. 

We find diffusion constants of Ki = 0.051±0.012 in the smaller box, and Ki = 0.073±0.026 
in the larger box, which agree within the statistical errors. Baryon number violation is about 
30 times less efficient than it would be if the diameter of the symmetry restored region were 
1 / awT, where the usual symmetric phase rate would be recovered . The sphaleron rate 



on the string is orders of magnitude larger than the rate in the broken phase, so the string 
is strongly enhancing baryon number violation; however, the rate is not quite high enough 
to say that sphalerons are fully unsuppressed along the string. 

We also find that k' is strongly dependent on 0o, and hence on the temperature. Going 
to fo/T = 1.4, it falls to a value of k' = 0.028 ± 0.008. Since 0(T) ~ <p{T = 0)(1 - T 2 /T 2 ) 
except very close to the phase transition temperature, and since T c ~ 100 GeV for a realistic 
scalar self-coupling and top quark mass, changing (po/T from 1.0 to 1.4 corresponds to a 
decrease in temperature of approximately 8%. Thus k' is quite temperature sensitive at 
Q /T~1. 
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2.3 Fat strings 

Next we return to the first example of subsection |2.1| , and investigate the possibility that the 
resulting string is "fat," with the Higgs field VEV suppressed in a large region surrounding 
the string. We would like to find values of the quartic couplings and the singlet VEV So 
defining the potential (||) which give the largest rate of sphaleron processes inside the largest 
possible strings. This requires that the doublet VEV <p(r) be close to zero near r = 0, the 
center of the string, and throughout a region of space large enough to allow unsuppressed 
baryon number violation. 

Before examining this model's suitability for baryogenesis, it may be worth pointing out 
that it cannot be obtained from supersymmetry, because an F-term contribution to the 
potential would give a positive IS 1 ! 2 !^ 2 coupling, whereas we need it to be negative. A 
D-term would give a negative coupling if the $ and S fields had opposite charges under a 
£7(1) gauge symmetry, but this situation is not compatible with having a fat string, for the 
following reasons. Certainly we should not identify the U(l) with weak hypercharge, since 
the S field would break it at a high scale. On the other hand, if the standard model Higgs 
field $ transforms under a new U(l)', the mass mz> of the corresponding Z' boson must be 
greater than 800 GeV, as discussed in section 2.1; however the string width is of order m^}, 
and so the string will be narrow, not fat, such as we wish to construct here. 

Also, if the U(l)' symmetry associated with S — > Se td is gauged, then for fixed So and 
X s , the resulting string is always thinner than if the symmetry is global. This is particularly 
true if the new gauge coupling g' satisfies g' 2 ^> \ s , as it must if A s is small and g' is to unify 
with the standard model gauge couplings at some high scale. As we will see, the case of 
small X s is the most promising. We also observe that the global string's fields approach their 
asymptotic values at r = oo slower than those of a gauged string. The scalar condensates 
have power-law behavior in global strings, 0(r) ~ 0(oo) — ar~ b , but exponential behavior 
in local strings, <p(r) ~ 4>{oo) — ae~ br . Therefore global strings also have a better chance, 
qualitatively, of having a large region of symmetry restoration. For these reasons, we will 
consider only global strings in the present section. 

While searching for the parameters favorable to baryogenesis, we must keep in mind a 
number of constraints: 

1. The vacuum must be stable, implying that 

4AA S > 7 2 . (14) 

2. The mass of the lightest Higgs boson must be above the experimental limit. While this 
limit depends on the relative admixtures of the $ and S fields in the light eigenstate, 
it is safe to say that m h > 75 — 80 GeV. Solving for the mass eigenvalues gives the 
constraint 

K = rn 2 h /$l 

= tyoz (x/y + y /x - y] (x/y - y /x) 2 + > 0.2 , (15) 
where we have defined the useful combinations 

x = \J\/X S 
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y = (S'o/$o)|t=t c ; Vo = OSb/$o)|r=o 



z 




1 = l/z, (16) 

and $ = 246/^ GeV = 174 GeV at T = 0, but it has a smaller value, v^o > T ~ 
100 GeV at the temperature where the bulk sphaleron rate first becomes negligible. 

Really the condition should be stronger still; we should also demand that either the 
lighter vacuum mass eigenstate must be almost purely $, or both eigenstates must be 
fairly light. Otherwise electroweak radiative corrections could be generated that would 
contradict precision tests. (Another way of evading precision electroweak test bounds 
would be to make the heavier mass eigenstate almost pure Higgs field, with a mass of 
~ 100 GeV; but this turns out to require parameters that do not lead to electroweak 
symmetry restoration in the string core.) However, in what follows we will be able to 
rule out string-mediated baryogenesis without bothering to enforce this condition. 

3. We require couplings to be perturbatively small, as otherwise the theory receives large 
radiative corrections, and actually cannot be defined without a UV regulator close to 
the energy scale of interest. Our condition (perhaps too generous) is 

A, \„ 7 < 1 • (17) 



4. Finally, a less obvious but crucial requirement is that the U(l)' symmetry broken by 
the VEV of S be restored at high temperatures. Otherwise the phase transition that 
should have given rise to the strings will never have occurred. This means that the 
thermal correction to the S field mass must be positive: 5m 2 = (2A S — 7)T 2 /6 > 0, 
and hence 

7 < 2A S . (18) 

This condition might be evaded by adding extra particles to the theory which would 
enhance the thermal mass of the S, about which we will say more below. 

The search of parameter space can be simplified by noticing that a rescaling of the string's 
radial variable r by r — > r/z 1 / 2 transforms the Hamiltonian density according to 

H(r; A, A s , 7 ) -> zH{rjz x l 2 , x, l/x, 7 ). (19) 

This implies that the solutions for the fields around a string can be obtained from the case 
where z = 1 by simply stretching the spatial size by the factor \j z 1 ! 2 . The other important 
measure of symmetry restoration, the smallness of h(r) at r = 0, is unaffected by this 
transformation. So it is sufficient to vary just the three parameters x, y and 7 of eq. ([H]) to 
generate all possible solutions. 

The most favorable case for baryogenesis occurs when 7 is pushed close to its upper 
limit of 2, coming from vacuum stability. Naturally one needs for 7 to be large since in the 
opposite limit of 7 = there is no coupling between the S and $ fields, hence no possibility 
for electroweak symmetry restoration in the strings. Just at the vacuum stability limit the 
potential develops a flat direction, hence a massless Higgs boson. The experimental limit on 
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the Higgs boson mass, together with Eq. ( ]T7|) , prevent one from taking 7 higher than about 
1.8. 

To make 0(0) small, one also needs to have So 3> $0 or y ^> 1. Otherwise the symmetry 
is only weakly or not at all restored in the core of the string. On the other hand, we need a 
wide string core to avoid the thin string case, which already proved unsuccessful in subsection 



27T| . This requires that X s Sq not be too large, hence A s must be small. To summarize, the 



most favorable situation for baryogenesis is when 

z~y>l; 7<2. (20) 

The field profiles in such a case which is favorable for baryogenesis are shown in figure 
|2|. There we graph $(r) and S(r) versus r for the couplings A = 0.35, X s = 8.8 x 10 -5 , 
7 = 9.4 x 10~ 3 , and S = 630 o (corresponding to x — 31, y — 63, z = 5.5 x 10 3 and 
7 = 1.7). The overall scale of the couplings is chosen to saturate the experimental bound on 
the Higgs mass. One notices that the region of symmetry restoration is very wide compared 
to the weak scale $0 1 = 10~ 3 fm. In this example, the half-width ri/2, where the Higgs field 
reaches a value of ($(0) + $o)/2, is nearly 20 times $0 • To get such a wide $(r) profile, it 
is essential that we use global rather than gauged strings. No radial field excitation in the 
broken phase is light ,[] and if the S field approached its asymptotic value exponentially, the 
region of symmetry restoration could not be wider than the largest inverse mass. 

Figure shows how $(0), 77/2 and m\ depend on x and y for the choices 7 = 1.7 (same 
as in figure 2) and z = 1. The crucial observations are that $(0) tends rapidly toward zero 
(maximum symmetry restoration in the string) as Sq/^q gets large, and the Higgs mass 

1 There is, however, an axial excitation-the axion-which creates cosmological problems for these models 
when the symmetry-breaking scale is low. 




Figure 2: Doublet and singlet Higgs field profiles around a string, for the parameters 
A = 0.35, A s = 8.8 x 10" 5 , 7 = 9.4 x 10" 3 , and S = 630 o , where $ = 174 GeV. These 
parameters do not allow thermal symmetry restoration of the S field. 
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increases most rapidly along the direction of increasing x ~ y, while the half-width varies 
relatively slowly. To increase 7*1/2 dramatically, as in figure one should go to large values 
of m\ and x ~ y, and then rescale all the couplings by z = 0.2 /ml, so as to decrease the 
Higgs mass to the point of saturating the experimental bound. This increases 7*1/2 by the 
factor z^ 1 / 2 without changing $(0). 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 

log x log x log x 



Figure 3: (a,b) Dependence of the Higgs field central value and half-width, 0(0), 7*i/ 2 
respectively, on the potential parameters x = \J\/X S and y = S /& , for 7 = 1.7 and z — 1. 
The logarithm is base 10. (c) Contours of the dimensionless Higgs mass m\ = m\/<&\ for 
the same parameters. 

However we have not yet taken into account the symmetry restoration bound, Eq. (TT8|), 
which does not allow one to make A s much smaller than 7. If we saturate this bound, we find 
that 7 = 2/x, forcing us to give up at least one of the two baryogenesis-favoring conditions, 
either that 7 ~ 2 or 1 > 1. If we give up the latter condition, so that x ~ 1 and logx ~ 0, 
figure 3 shows that m^ does not grow with y, giving no leverage to rescale z and hence 
to widen the Higgs field profile. Indeed, eq. (|T5| ) shows that in the fixed-x, large-?/ limit, 
ml = xz(2 — 7 2 ). Thus the smallest we can make z is 0.2/(2 — 7 2 ), which is of order unity 
or greater if 7 is close to 2. The region of symmetry restoration is of order the weak scale in 
this case, which is too small to comfortably contain a sphaleron. In fact, the large-?/, x ~ 1 
case gives the thin strings we have already considered above. And if we choose x ~ 1 and 
y ~ 1, then $(0) is not close to zero; the Higgs field does not lose its condensate even in the 
core of the string. 

If, on the other hand, we try to keep x large while decreasing 7, the electroweak symmetry 
restoring effects of the string rapidly diminish, as shown in figure 4, which is the same as 
figure 3 but for the smaller coupling 7 = 0.1. As the figure shows, it takes a very large value 
of y (So) to suppress <3>(0) in this case; however ri/ 2 decreases with y for fixed x, while ml 
reaches an asymptotic value. Taking for instance x = 20, the largest value allowed by the 
symmetry restoration bound at this value of 7, we get a string of radius 77/2 < 3.5/0 o and 
4>(0) = O.10o when y = 100. Both quantities decrease monotonically with y, which must 
be large in order to suppress the Higgs field in the string core, but the radius over which it 
is suppressed shrinks at large y. For this value, x — 20, y = 100, and 7 = 2A S , we find a 
sphaleron energy 0.492 times the broken phase value. The string cuts the sphaleron energy 
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Figure 4: Same as figure 3 but with 7 = 0.1 



in half, so the sphaleron rate will be substantially faster than the broken phase rate, but still 
much slower than the symmetric phase rate, because of the Boltzmann suppression factor 
e -s sph /T_ rj^ e situation does not improve if we vary y in either direction. 

It is possible to contrive to evade the symmetry restoration bound, eq. fll%|). One could 
add another scalar field / with a positive quartic coupling r y'\f\ 2 S*S, whose sole purpose is 
to increase the thermal mass of the S field. The constraint is then relaxed to 7 < 2A S + 7' 
if / is complex and 7 < 2A S + Y/2 if / is real. But adding / introduces new problems: 
A s receives radiative corrections of order ('~f') 2 , which requires fine-tuning of the ultraviolet 
couplings to preserve the hierarchy X s <C |7|. Furthermore we must either tune mj < 1 
KeV, or give the / particles a direct decay channel by adding a cubic term of form /|$| 2 
and making rxij > 4m;,; otherwise they leave an excess relic abundance which overdoses the 
universe. 

It appears possible, then, to allow for "fat" strings with large regions of symmetry restora- 
tion, but it requires a highly-tuned hierarchy of couplings in a very artificial model. Aside 
from this possibility, "fat" global strings do not appear to provide an effective means of 
restoring electroweak symmetry in a suitably large string core to allow efficient baryon num- 
ber violation. 

To conclude this section, there are two nonsuperconducting string models which allow 
significant baryon number violation along strings. The case where an extra U(l)' is broken 
by a complex scalar works if the SM Higgs field has a charge under the U(l)' which is 
incommensurate with the charge q s of the scalar responsible for breaking the U(l)' symmetry. 
The baryon number violation along such a string is much faster than in the broken phase, 
but is not completely unsuppressed. It is also possible to violate baryon number efficiently 
in a global string model if the coupling between scalar and Higgs is negative and if there is a 
large hierarchy of scales between the Higgs self-coupling, the Higgs-singlet coupling, and the 
singlet self-coupling. However it is also necessary to add a real scalar with cubic couplings 
and to tune them in a way which is unstable to radiative corrections; thus the model appears 
to be rather unnatural. 
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3 Generation of CP asymmetry inside the string 



As we saw at the beginning of the last section, the baryon number production is proportional 
to the product of the sphaleron rate along the string and a CP-odd sum of particle chemical 
potentials, averaged over the volume in which the sphalerons occur along the string. So 
besides the sphaleron rate, discussed in the last section, we also need to know how the string 
biases particle populations, so as to create an asymmetry between particles and their CP 
conjugates inside the string. 

Chemical equilibrium together with CPT ensures that CP-odd combinations of chemical 
potentials must vanish, so the string only generates such particle distributions if it drives 
the plasma out of equilibrium, which occurs if the string is moving relative to the plasma. 
The main goal of this section is to show that, if the string moves slowly, any CP violating 
chemical potential it generates depends quadratically on the string velocity v s . The fact that 
it goes as v 2 s and not v s is easy to understand; a string has no natural "front" and "back" 
sides, so there is no distinction between it moving forward or backward; hence the chemical 
potentials averaged over the string volume must be even functions of v s . Assuming only that 
the chemical potential on the string has analytic behavior in v s , we get that // ~ v 2 for small 
v s . It remains to check whether this picture is really right, and to see how big the coefficient 
is. 

We will treat this problem in three different ways. First, we will consider an unrealistic 
model where all parameters are pushed to the extreme values that give the most favorable 
outcome for baryogenesis. This will be used in the following sections to make a robust 
argument against the viability of string-mediated baryogenesis. We will then make more 
realistic estimates of the CP-violating chemical potentials, ficp, using a purely quantum 
mechanical treatment for the particle-string interaction, as well as a semiclassical one. It 
will be seen that both of these give smaller predictions for the size of the CP violating effects 
than does the toy model. A crucial aspect of these results is the dependence of ficp on v s , 
the string velocity. In particular, for small v s we will show that Hep oc v 2 s . 

3.1 Model-independent upper limit on /icp 

CP violation at the string wall comes about when the probability for particles to reflect 
from the wall differs from that of antiparticles. The most extreme situation imaginable for 
maximizing this effect would occur if the string walls totally reflect particles, and perfectly 
transmit the corresponding antiparticles. We will make this radical assumption now, in order 
to get a robust upper limit on /i<7P> hereafter simply denoted by /j,. Let us now consider 
how this affects the particle distributions inside the string, where baryon number violation 
can occur. Since the antiparticles do not see the wall, they are continually refreshed by the 
equilibrated population outside of the string, so their distribution functions remain thermal. 
We thus have 

f~ = e E/T ± I = ff ' ( 21 ) 

where // denotes the equilibrium Fermi population function. 

The particles, on the other hand, have in addition the effects of the reflections, which in 
the rest frame of the plasma causes their distribution functions to be boosted by the string 
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velocity (taken to be in the x direction). It is also possible that they acquire a chemical 
potential, but other departures from equilibrium are damped by collisions. Assuming only 
that collisions with plasma particles are less frequent than collisions with the walls of the 
symmetry restored region, which is reasonable for the relatively narrow strings considered 
in this paper, we get 

/+ = e ^-J + » )/T±v ^ = (1 - ^ ■ (22) 

For definiteness we will henceforth assume the particles are fermions, and denote them by 
F, since only fermions enter in the anomaly equation and bias baryon number violating 
processes. 

The generation of the chemical potential [i is due to annihilation [FF — > GG) and pair- 
production processes, (GG — > FF), involving gauge bosons G, for example. The size of the 
equilibrium value of fi which is generated can be deduced from the Boltzmann equations for 
the F and F distributions: it is that value of ji which causes the production rate of particles 
via GG — > FF and the destruction rate from FF — > GG to cancel each other. This is the 
right condition because particles inside the string can never get out; if the production and 
destruction rates did not balance, the particle number would adjust until it did. Let /& = 
[e E l T — 1) _1 denote the thermal distribution functions for the gauge bosons, M. the matrix 
element for annihilation or pair production, and dli = 5(p\ +p 2 ~P3 ~Pi) Y\t=i d 3 Pi/(167i 3 Ei) 
the phase space measure. The condition that the collision term vanishes is 

J dU\M\ 2 [/ + ( Pl )/_(p 2 )(l + f b (p 3 ))(l + f b ( P4 )) 

-/6(ps)A(ps)(l - /+(pi))(l " f-(P2))] = 0. (23) 

Expanding f± in powers of /i and v s , and using that E-y + E 2 = E 3 + E4, we find after menial 
algebra that the term in brackets involving population functions is 

[//(Pi)//(p 2 )(l + fM)(l + hipi))] + v s V -f - ||J + |^(l - 2/^Pi))) , (24) 

plus corrections of order /i 2 , v 3 ., or fiv s . Since the integration measure contains an average 
over angles, p z averages to zero and pi averages to (l/3)p 2 . It is because p z averages to zero 
that no v s term appears in the expression for ji. It is also clear that for all higher order 
processes, the term proportional to v s will involve an integration over a vector quantity and 
will vanish on angular integration; so the v 2 dependence of fi is not an accident of expanding 
to leading parametric order. 
Now, defining 

(x) IdU\M\ 2 f f ( Pl )f f (p 2 )(l + / 6 (p3))(l + f b (p,))X 

{ )M ~ /dn|A<|V/(Pi)//(P2)(i + A(ps))(i + /6(P4)) ' 1 ] 

then the result for the chemical potential can be simply expressed as 



v 



2 



l l^^l T - V 



1 ( ^-3tF^Tt) I (26) 
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Notice that the result is independent of the overall strength of the interaction, Ai. This 
will be true so long as the interaction is faster than the Hubble expansion rate, so that it 
remains in equilibrium. The answer does involve the energy dependence of A4, but the fact 
that fi oc v 2 does not. 

To give a concrete estimate for the size of /x, we can simplify the integral by considering the 
leading logarithmic behavior in the gauge coupling constant g, which comes from t-channel 
exchange at low momentum transfer, where M. ~ l/t. Using the same approximations as in 



Appendix A of ref. | 27fl , we obtain 



» = j^v 2 s T 0.12v 2 s T, v s <l, (27) 

in the small- velocity limit. The value of /i, relevant to baryogenesis, is the sum of 1/2 this 
quantity, over all left handed doublets coupled to SU(2) which scatter from the wall in this 
way. 

While we have only performed the calculation in the small v s case, where it is valid to 
make an expansion about equilibrium population functions, we expect that, for relativistic 
strings, 

H<T, -y v v a ~l. (28) 

since parametrically there is nothing that can enhance /i to the point of getting a highly 
degenerate gas of CP-asymmetric particles. In fact, it will be seen that the arguments to be 
made below would not be invalidated even if \i diverged like j v in the ultrarelativistic limit. 



Eqs. (|27|-p8|) are sufficient for making our estimate of the maximum string-mediated 
baryon asymmetry possible, as will be done below in section 4. To bolster and complete the 
argument, however, we shall give some more realistic computations of the CP asymmetry in 
the following subsections. 

3.2 Initial chiral flux in a realistic model 

We will now consider the generation of a CP-violating particle flux in a specific theory. The 
simplest example is a two-Higgs-doublet model including |28| a complex mass term m?${$ 2 



and quartic interaction /i($|$2) 2 - The unremovable CP-violating phase 9q is defined by 
m^/i 1//2 = \m^h}l 2 \e ld ° . In the string wall, the two VEVs can be parameterized as $i(r) = 
0i(r)e l6 '( r )/ 2 and — (f)2(r)e~ l9 ^ r ^ 2 , where the <pi functions are real and 59 = 9(oo) —9(0) 

is of order 9q. Solutions for 9(r) have been found for domain (bubble) wall backgrounds in 



ref. [28], where 4>i(r) has the form tanh(r/r ), and we expect the solutions to look similar in 



the case of strings since the shape of our Higgs doublet profiles are qualitatively the same. 



We will therefore borrow results from ref. [28| to make the following estimates. This is done 



to be as quantitative as possible, but our main conclusions will not depend on fine details of 
this particular model, but rather on parametric dependences that we expect to hold in any 
model. 

If each fermion couples only to one of the two Higgs doublets, as is beneficial for avoiding 
flavor-changing neutral currents, then they will have CP-violating reflections from the string 
wall as a result of the spatially varying phase 9(r). Those particles with the largest Yukawa 
couplings will be reflected most strongly. For simplicity we will model the string as having 
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a square rather than round cross section, moving through the plasma perpendicularly to 
one of the sides (□ — > x) with velocity v s = v s x. While this approximation is far from 
perfect, it can only change the answer by a geometrical factor of order 1. It allows the use of 
planar wall calculations for the difference in reflection probabilities, which can be adequately 
parameterized in terms of the fermion momentum perpendicular to the string wall (which 
we will take to be p x ) by an exponential, 

AR(p x ) = A56e~ Px/Sp (29) 

(the linear dependence on 56 being correct for 56 <C 1), where A and 5 P depend on the mass 
mj of the fermion outside the string and the width of the Higgs field profile. The dependence 
of 5 P was found to be 



5 P _ 
m f 



—1.1 \n(mfri/ 2 ) — 0.54, mfiji < 0.3; 

0.19 (mjr^)- 1 - 2 , 0.3 < m f r 1/2 < 0.7; (30) 

0.15 (m/r^ 1 ' 8 , m f r 1/2 > 0.7; 

It should be noticed that, although large values of T\j 2 are desirable for enhancing the rate of 
sphaleron interactions in the string, they at the same time suppress the quantum reflections 
of heavy particles. For this reason lighter particles like the r lepton can be more important 
than top quarks despite their smaller Yukawa couplings. 

The amplitude of the asymmetry, A, depends not only on mfTi/ 2 but also on w?^jvp? h . 
For m/ri/2 ~ 1, A is also of order 0.1, but for larger m/ri/ 2 it decreases exponentially, 
A = e~ cmVl / 2 , where c ranges from 0.5 to 2 as m 2 ^/m\ goes from 0.3 to 8. 

Once AR is known, the next step is to compute the flux of CP asymmetry going into the 
string. The flux of CP current per unit length of the string is given by 



J 



1 r d 2 p ± dp x 

2 J (2vr) 3 E 



AR(p x )(f( lv (E + v s ^jpl-m))/T) - f{ lv {E - v aPx )/T)), (31) 



where / is the distribution function of the particles, boosted by ±v (the wall velocity) 
depending on whether they are approaching from the right or the left, and 7„ = (1 — f 2 ) -1 / 2 . 
Making the Maxwell-Boltzmann approximation for the distribution functions, and ignoring 
the dependence on mj (since the particles are assumed to be relativistic) we can easily 
evaluate J to be 

J = v a A66% 1+ W T 5 (32) 

vr 2 (l + 2 lv 5 p /T+(5 p /T) 2 ) 2 1 J 

In another theory the form of J might differ, but it must always be proportional to the string 
velocity, v s , because at v s — everything is in equilibrium. Furthermore, it is bounded by 
J ^ VsT 3 /^, which is the forward-backwards difference in total incident particle numbers. 
This bound would be saturated only if the string were a perfect reflector of one CP state and 
a perfect transmitter of the other; any realistic wall will give currents that are much smaller. 



3.3 Subsequent diffusion of the source 

Now we are ready to use the current as a source for the diffusion equation to find the steady- 
state density of chiral asymmetry at the center of the string. It must be understood that the 
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flux J actually represents a positive flux entering the front wall and going into the string, 
and an equal and opposite flux leaving the string through the front wall and entering the 
plasma. To apply the diffusion equation, it must be assumed that these fluxes penetrate the 
plasma by some distance A (the mean free path) before diffusive behavior sets in. Otherwise 
the diffusion equation sees only two equal and opposite, hence canceling, local densities at 



the string wall, giving no net source of chiral charge. As explained in references p8| , |29| , this 
can be modeled by a source term in the diffusion equation which is a difference of two delta 
functions: 

dn ( d 2 d 2 \ 

fr- D {d^ + w) n + Tn = J {5{X + A ~ Vst) - 6 ( X ~ X - «'*)) ( 33 ) 

Here y is the transverse direction along the cross-section of the string, D is the fermion's 
diffusion constant, and T is the rate of chirality-damping processes, like spin- flip interactions. 

However, the above equation gives only the contribution from the front wall. There 
is another contribution from the back wall which has the opposite sign. Taking this into 
account, and changing coordinates to the rest frame of the string, where x — y — denotes 
the center, we find that the solution for n can be expressed as a sum of four contributions, 
two for each wall: 

n{x,y) = j[N(x + w/2 + X,y)-N(x + w/2-X,y) 

+N{x - w/2 + X, y) - N(x - w/2 - A, y)] , (34) 

where w = 2r\/2 is the width of the string and N is the Green's function of the 2-D diffusion 
equation. The first line represents the front wall contribution and the second is that of 
the back wall. Although the explicit signs in the equation are the same for both walls, the 
contributions to the chemical potential at the center of the string are of opposite sign; it is 
the second and third contributions in the equation which are injected inside of the string. 
It is only because the string's motion leads to an asymmetric solution N of the diffusion 
equation that the expression for n is nonzero at the center of the string; in fact as the string 
velocity v s goes to zero N becomes a symmetric function and the average over the string 
interior of the term in brackets goes to zero linearly in v s ; this is whence the second power 
of v s arises. 

Using the Green's function for the two-dimensional diffusion equation, the individual 
contributions are given by 

N(x, y) = — r — f ^ ^' e -(«/ 4 ^')((^') 2 +(,-/) 2 )-rx7^ (35) 
AnD Jo x' J-w/2 

which satisfies a two-dimensional diffusion equation like (|33D with source term 

/■to/2 

5{x) dy'S(y-y'). (36) 

J-w/2 

The derivation of eq. ([35]) is not completely obvious. For a stationary source, the Green's 
function for the 2-D diffusion equation is G(x,y,t) = t _1 exp(— (x 2 + y 2 )/4Dt — Ft). Now 
we are considering a moving source, whose trajectory can be described by x' = vf. At 
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a given position (x, y) and time t, an observer sees a contribution from when the source 
was at position x' at the earlier time t', proportional to / dy'G(x — x',y — y',t — t'). We 
must integrate all these contributions from t' = — oo to the present time t: N(x,y,t) oc 
ftoo dt' J dy'G(x — x',y — y', t — t'). Here N is still in the rest frame of the plasma; replacing 
x —>■ x + v t takes us to the rest frame of the source, so that x now measures how far one is 
in front of the source, and the solution becomes time- independent. The change of variables 
x' = v(t — t') in the integral gives eq. (|35|). As for the normalization, it is chosen to insure 
that the rate of particle creation by the source equals that of particle loss by decay when 
the system reaches a steady state: 

r J J N(x,y) dx dy = Jw. (37) 

At the center of the string (y = 0), where the effect is maximized, the answer simplifies 
somewhat, and can be expressed as 

—vx/2D J 

N(x) = N(x, 0) = =■ / ^= etf(w/V^)e- x2/z -( v/4D) ( v/4D+r/v)z . (38) 

ADy/TC JO V Z 

Substituting this result into eq. (0) and evaluating it as a function of x, we find that the 
chiral density in the string is almost antisymmetric about the core, but is skewed slightly by 
the nonzero velocity, which distinguishes front from back. This behavior (and the fact that 
we have idealized the string by four delta function sources, eq. (^)) is shown in figure [5[ 
Because n(x) varies significantly inside the core of the string, and the rate of baryon number 
violation at position x is proportional to n(x), we should average it over the core, i.e., the 
region x G [—w/2,w/2\: 

1 W2 

n — — I n(x, 0) dx. (39) 

W J-w/2 

We should also average over y, but n is a smooth function of y, peaking at y — 0, so using 
y = only slightly overestimates the average. This approximation, in any case, errs in favor 
of the mechanism we are ruling out. 

In Figure ^ we show how the spatially averaged chiral asymmetry at the center of the 
string varies with the string velocity, for several choices of parameters. Set 1: D — 6/T, 
A = D, T = T/100, w = 6/T, 5 P = 0.1T, which is a realistic choice for quarks [[30(1 . Set 2: 



same as set 1 but with D = X = 100/T, which might be reasonable for leptons. Set 3: same 
as set 1 but with 5 P = T, to show what happens when the bubble surface efficiently reflects 
particles with higher momenta. Set 4: same as set 1 but with w = 100/T, to show what 
happens for a very thick string. In all examples we have assumed that T = 100 GeV, which 
is the temperature where the electroweak phase transition is known to occur, and that CP 
violation is nearly maximal, 59 — 1. 

In every case the number density, and hence the chemical potential (related by n = fiT 2 /6 
times the number of spin degrees of freedom), varies quadratically with the string velocity, 
/i/T = KVg, with K significantly less than 1. One power of v s arises from the injected flux, 
while a second power comes from a front-back cancellation. 

To understand why the coefficient of the v 2 s law is so small, it is important to realize 
that the string core width is typically not much larger than the quark diffusion constant, 
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Figure 5: The chiral asymmetry as a function of distance from the center of the string, for 
the parameters D = 6/T, A = D, T = T/100, w = 6/T, 5 P = 0.1T, and v = 0.1. 

and so an asymmetry generated on the string fairly rapidly diffuses away from it into the 
plasma at large. Since baryon number violation is only efficient right on the string, any 
chiral asymmetry outside the core is wasted, as far as baryogenesis is concerned, which helps 
to explain why the coefficient of \i is so small. Even if we took the flux to be maximal the 
coefficient K would still be of order 10~ 2 . This is to be contrasted with the situation in the 
conventional scenario with a first order phase transition; here the wall is planar, so there is 
only one direction to diffuse away from it, and baryon number is efficiently violated in one 
of the half-spaces. 

It should be mentioned that the f s -dependence of the chemical potential would be 
different for a very thick string. When the string width satisfies both r > D/v s and 
r > T 3 /r sp h ~ 10 6 /T (where r sp h ~ 10~ 6 T 4 is the rate per unit volume of sphaleron 
interactions in the symmetric phase) then the baryon number asymmetry from the front 
side vanishes; any baryons produced fall into the string and are inevitably destroyed by 



sphalerons |T(|. In this case the baryon number production can scale linearly with v s (al- 
though at small enough v s , the r > D/v s condition will eventually break down). However, 
the nonsuperconducting strings we have considered have radii nowhere near this large, and 
/iocifj holds up to v s ~ 1. 

Our treatment becomes less good for v s approaching 1, but we do not expect a large CP 
asymmetry in this case either. At such large velocities, the flux impinging on the back side 
of the string can be ignored compared to that which hits the string's front side and escapes 
through the back. In the string rest frame, this incident flux is composed of very hard 
particles. The reflection probability for hard particles drops very quickly with momentum, 
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Figure 6: The chiral asymmetry averaged over the string cross section, as a function of the 
string velocity, treating particle reflections from the string quantum mechanically. See text 
for description of parameter sets. 



as was mentioned above, and so therefore does the CP-violating difference in reflection 
probabilities. In any case, even in a scaling string network where friction is negligible, almost 
none of the strings are traveling ultrarelativistically; the mean value of v 2 s (with the average 



weighted by the energy and not the length of string) is less than 1/2 [32]. For a scaling- regime 
network what is relevant is an appropriate average of fi(v s ) over the velocities represented by 
the scaling network, for which we expect fi ~ T with a constant of proportionality roughly 
of the same order as K the coefficient in the v\ law discussed above. 

3.4 WKB estimate of the chiral asymmetry 

It has been pointed out that when defect boundaries are sufficiently thick compared to the 
mean free path for particles to scatter in the plasma, a semiclassical treatment of chiral 
asymmetry generation is more appropriate than the quantum mechanical one given above 



3T| , and gives a larger answer. The distinction arises because CP violation in the defect 
creates a classical, CP-violating force on particles as they cross the wall. In addition, the 
scatterings themselves can be CP-violating, giving rise to what has been called "spontaneous 
baryogenesis" in earlier literature. 

If d{x) is the background CP- violating phase, such as appears in the two-Higgs-doublet 
model discussed above, then these two classical effects generically give rise to a source term 
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Figure 7: The chiral asymmetry at the center of the string as a function of v s , using the 
WKB approximation, and the parameters D = 6/T, Thf = T/100, w = 6/T, 89 = 1, yf = 1. 



where is the Yukawa coupling of the fermion to the Higgs field. The primes denote spatial 
derivatives, and T^f is the rate of helicity-flipping interactions. The first term is due to 
the classical force, and the second is the origin of spontaneous baryogenesis. The size of 
the chemical potential thereby generated can be roughly estimated just as in the previous 
section if we approximate 9 as a step function, so that 9' = S9S(x). Taking into account 
the separate contributions from the front and back walls of the string is accomplished by 
replacing the above source term by 



The minus sign occurs because the shape of the 9(x) profile on the back wall is the mirror 
image of that on the front wall, and S(x) involves an odd number of derivatives of 9(x). This 
behavior is in contrast to the quantum mechanical calculation, where front and back walls 
gave same-sign contributions. This difference between the quantum and classical treatments 
stems from the fact that quantum reflection probabilities are the same for particles crossing 
a potential step, regardless of the direction of motion, whereas the sign of the classical force 
exerted on the particle by such a step does depend on what direction it is going. 




(40) 



S 2 ( x ) = S(x + w/2) - S(x - w/2) 



(41) 
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With the delta-function approximation for 6', we can immediately write the spatially av- 
eraged solution for the chiral asymmetry in the string, in analogy to the quantum mechanical 
treatment of the preceding section, n is given by eq. (|39|), using 

<x) = v s 69 y -lf (DN"(x>) + T hf N{x'))Cj:_ w j; , (42) 



and the definition of N(x) given by eq. (P8|). The results for the previously-described param- 
eter set 1 (D = 6/T, Thf = T/100, w = 6/T, 59 = 1), assuming yf = 1, are shown in figure 
[7|. The two terms, due to the spontaneous baryogenesis and classical force effects, respec- 
tively, are shown separately. The classical force contribution is larger than the spontaneous 
baryogenesis contribution, and is also significantly larger than the quantum mechanical pre- 
diction for the same set of parameters. The v 2 dependence is evident over most of the range 
of velocities. To convert from density to chemical potential, one uses 



for a fermion with g helicity states. 



4 Efficiency of baryogenesis by a string network 

We have now seen how efficiently baryon number is violated along cosmic strings. Namely, 
the production of baryons per unit length of string is 

^ = 1.5^T 2 V^, (44) 
dLdt w Y T 

where the \ii are the chemical potentials for the ith species of left-handed (SU(2) doublet) 
fermions. Henceforth we will write YaM simply as fi. We have already seen that Ki < 0.1, 
falling fairly rapidly as the temperature drops, and that \i varies as v 2 , the square of the 
string velocity, when v s is small. To estimate the baryon number produced, we still need 
to know the density and velocity of the string network, however. These quantities depend 
crucially on one model-dependent parameter of the strings, namely their tension. 

Let us denote the string tension by r, the typical string's radius of curvature at the 
electroweak epoch by R, and the typical velocity of propagation of a string, in the rest frame 
of the plasma, by v. On dimensional grounds we have that 

v~HR, (45) 

with H the Hubble constant at the electroweak epoch. 

The string experiences a velocity-dependent friction as it moves through the plasma, 
arising from particles bouncing off the electroweak-symmetry-restored core of the string. 
For small velocities one can define a frictional constant rj, such that the force per unit length 
of string is dF/dL = —rjv. If w is the width of the symmetry-restored region, then 

ri-wnb-J-g^T, (46) 
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where r/b is the friction constant of a moving electroweak bubble wall, F/area = — vrfo, and 
its parametric dependence, rj^ ~ g 3 0Q_T, was argued in ref. If we take O ~ T, as 

needed to prevent baryon number violation in the bulk from erasing the baryons produced 
by strings, then rj ~ g 2 T 3 ; if instead we make the parametric estimate appropriate right 
after a weak phase transition, 0o ~ gT, we get rj ~ a^T 3 . This is the friction on the string 
just from its dragging a region of restored electroweak symmetry through the plasma; there 
may also be additional sources of friction, so the r] given above is a lower bound. 

The force (per unit length) pulling a string forward is given by t/R. If v <C 1 then inertia 
is not important and we can equate the forces, and using eq. (|5|), obtain the estimates 

i=„ - V ~^L,R~^L-, (47) 

valid whenever the resulting velocity is small. This situation is called the friction- dominated 
regime. On the other hand, if the above estimate gives v > 1, then friction is not important; 
instead the string evolution will be dominated by the tension, expansion of the universe, and 
loss of closed loops due to self-intersections. In this scaling regime, R ^ \/H and v ~ 1, 
with a mean squared value (v 2 ) < 0.5 |32|]. The crossover tension between the two cases for 
T ~ 100 GeV is r > (10 9 GeV) 2 . 

To obtain the rate of baryon number production per unit volume from eq. (fHj), we must 
calculate the total length of string per unit volume. One expects on dimensional grounds 
that the mean separation between strings is O(R), so the length of string per unit volume is 
l/i? 2 , and the baryon density production rate is 

dN B / 1 \ 2^2^ (a q\ 

^ t ~\W Kia ™ T T- (48) 

As was discussed in the last section, at small velocities the chemical potential behaves like 
/i = v 2 KT, with K a pure number depending on <fio/T but not on v. It was shown that K 
cannot be larger than 0(1), and realistically is much smaller. Similarly, in the large- velocity 
case, we expect p < T, say, /i = K'T. 

The two regimes of string evolution can now be considered in turn. First consider the 
friction-dominated case. Using fi = Kv 2 T and eq. (|45| ) in ( ^8|) gives 



dN B . .. .2 , • , • 

dVdt 



K\a w KH T . (49) 



This expression must be integrated from the time when the Higgs condensate first becomes 
large enough to preserve baryon number in the broken phase, 0o — T, to the present time. 
The right hand side depends on time through the temperature, d(\nT)/dt ~ H. The most 
strongly temperature-dependent term is k/, which decreases quickly with falling tempera- 
ture. Even making the generous assumption that it takes one Hubble time for K\K to fall 
sufficiently to cut off baryon number production, we still get a final baryon density of order 

Nb - ^Ka 2 ^ ~ KtKai— . (50) 



VT 3 1 W T 1 w m p i 

The last step follows from the Friedmann equation, H 2 = (8n/3)(p/m 2 ] ), using p ~ T 4 , so 
H ~ T 2 /m p \. In this expression, all dependence on the density of the string network, i.e., on 
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l/R, has dropped out. This happens because, although a denser network of strings produces 
baryons along a greater total length of string, these strings are moving more slowly and are 
therefore less efficient at baryon number production. 

The above argument neglects the effects of small closed loops of string, yet the string 
velocity in a friction-dominated network can become large for a loop just as it collapses. 
One may wonder whether such loops can enhance the baryon production. To answer this, 
consider a loop initially of radius r = R, the characteristic curvature length of the network, 
which subsequently shrinks. Its velocity is v(r) oc 1/r, and the radius evolves as dr/dt = —v. 
Hence r oc (to—t) 1 / 2 , with to the time of final collapse. The baryon production is proportional 
to t 

Lv 2 dt = 2tt f ° rv 2 dt oc f ° [t - t)~ 1/2 dt , (51) 



which is well-behaved at the upper limit of integration and dominated by r ~ R. Hence 
the collapse of loops contributes approximately as much to the baryon asymmetry as the 
evolution of the rest of the network, parametrically no more. 

Putting in the most optimistic values, k\ ~ 10" 1 , K ~ 10~ 2 , T/m p \ ~ 10 -17 , we 
find a baryon asymmetry about 13 orders of magnitude smaller than the physical value, 
N B /(VT 3 ) ~ Ng/Nj ~ lO" 10 . 

If the string network is in the scaling regime, the previous argument changes as follows. 
Rather than // = Kv 2 T, with K < 1, we have fi = K'T with K' < 1. Also, R ~ l/H 
rather than R ~ v/H. The final expression is the same with K' substituted for K, and the 
generated baryon number is still too small by at least 13 orders of magnitude. 

These estimates are optimistic, in the sense that Ki falls quite quickly with temperature, 
so there is much less than one Hubble time for baryon number to be created; moreover 
K is realistically much less than 10~ 2 . The treatment of the defect network is somewhat 
pessimistic, though, since numerical studies show that the mean separation of cosmic strings 



is closer to 1/10 of the Hubble length, rather than the full Hubble length [32fl. This could 
enhance our estimate of the baryon asymmetry by a factor of 10 2 . Also, the coefficients of 
order unity in the Friedmann equation imply that H ~ 10 T 2 /m p \, winning us an additional 
factor of 10. But this still leaves the mechanism too weak by 10 orders of magnitude, even 
using the most optimistic assumptions. 

Interpretation: Free Energy- 
There is another way of understanding our result, by considering the available free energy 
in a network. Think of the out of equilibrium requirement of baryogenesis as the statement 
that a baryogenesis mechanism must be a way of converting available free energy into baryon 
number. 

The available free energy density in the conventional electroweak baryogenesis scenario, 
with bubble walls converting supercooled symmetric phase into broken phase, is the free 
energy difference between the phases, which is (T c — T nuc )AS, the temperature drop between 
equilibrium and nucleation times the entropy difference of the phases. The numerical value 
is parametrically 0{a\ v T A ) for A ~ g 2 , but it is always at least 10~ 3 T 4 when the phase 
transition is strong enough to prevent subsequent erasure of the baryon number. 

For the string network, the available free energy is the free energy contained in the net- 
work, which is the length of string times the tension (times the mean value of the relativistic 
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7 factor for the scaling network, but this is order 1). For the friction dominated network, 
the density of the network is ~ 1/-R 2 , which together with Eq. (fI7D gives a free energy 
density ~ Hr] <T 4 (T/m p \), independent of the tension. Since the baryogenesis mechanism 
is the same or similar as that in the conventional first order case, it would be surprising 
if this much smaller available store of free energy were able to produce anything like the 
same number of baryons. For the scaling network the available free energy is larger, though 
observational constraints require it to be less than 10 _6 T 4 . However, most of the free energy 
loss is not due to friction against the plasma. In fact, the amount which is being lost to 
heating the plasma directly is independent of the tension of the network, since the width of 
the symmetry restored region does not depend on the string width. In fact the free energy 
loss to the plasma is about rj times the area swept out by the network in a Hubble time, 
~ Ht], which is parametrically the same as for the friction dominated network. In either 
case, the network is just not capable of generating anywhere near the amount of departure 
from equilibrium needed to produce enough baryons. 

The situation might be different for a superconducting network, but the above reasoning 
suggests that even for this case, baryogenesis can only get close to the efficiency of the con- 
ventional first order mechanism if the free energy of the network is close to the observational 
limit and much or most of the network's energy is lost to the plasma. 

5 Conclusion 

String networks do not provide a viable scenario for electroweak baryogenesis if the strings 
are not superconducting. While they are capable of generating baryons, they do so far too 
inefficiently to account for the existing abundance. 

The first problem is that baryon number violation along a string, while certainly faster 
than in the broken phase, is not entirely unsuppressed. The core where electroweak symmetry 
is restored is never wide enough to allow symmetric phase sphaleron-like events; at best it 
brings the sphaleron rate to 1/30 of what it would be in the unsuppressed case, i.e., if the 
core region was 1 / a w T in diameter. 

Given the actual width of the region of symmetry restoration around a cosmic string, 
we find that if the network is in the scaling regime, it violates baryon number in far too 
small a total volume for it to add up to the observed abundance. If the network evolution 
is friction-dominated, the network is denser, but the denser is the network, the more slowly 
moving the strings are. This in turn makes them very inefficient at creating baryon number, 
since the production rate goes as v 2 for small velocity. In either case the generated baryon 
number is at least 10 orders of magnitude smaller than the observed abundance. 

The case of superconducting cosmic strings merits more careful study. The main challenge 
here is finding a realistic estimate of the hypercharge current carried by a typical string at 
the electroweak epoch. We have not attempted to study this problem here. However, our 
results make it clear that unless the region of symmetry restoration proves to be much wider 
than l/g 2 T, the superconducting scenario will also fail. 
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A Cooling to a saddle point 

Here we discuss an algorithm for numerically finding a saddle point of a Hamiltonian on 
a many- dimensional space, which we used for determining the energy for a sphaleron in a 
string background in Section |2]. 

Suppose that we have a Hamiltonian system with degrees of freedom and Hamiltonian 
H(<&). We require only that the coordinates $ are continuous. They could for instance be 
real numbers, or members of a compact manifold. 

Ordinary gradient flow cooling of the fields means evolving the fields under a cooling time 
r according to 

_ dH(*(T)) 

dr ~ d$ a • {b ' 

When the $ take values on a manifold the derivatives are understood in terms of the tangent 
space at the current location. A typical discrete implementation of this algorithm is 

<Mr + A) - $ Q (r) = -A ^|<- T)) . (53) 

When $ is a point on a manifold, the right hand side is a point in the tangent space of 5>(t) 
and we determine $(r + A) by starting at 3>(t) and following the geodesic curve of starting 
direction and length indicated by the right hand side. 

Provided that we make A small enough, this algorithm will reduce the energy until it 
hits a local minimum.^ "Small enough" means that A < 2/cj^ ax , where is the largest 
eigenvalue of the matrix of second derivatives of H. In many cases we can actually compute, 
or at least bound, u^ ax over the entire space, and ensure stable dissipation. 

What if we want an algorithm which will flow to the nearest extremum, even if it is a 
saddle point and not a minimum? If we can guess a starting point which is fairly close to 
the desired extremum then it makes sense to analyze algorithm behavior in terms of small 
deviations about the extremum. Write the fields as $o + 8<& a , and expand the Hamiltonian 
to second order in 5Q a ] 

H~H{%)+H af3 6<S> a 8<S>f3. (54) 

The matrix H a p is diagonalized by the basis of eigenvectors £ Q , with eigenvalues X a . It is 
convenient to write the distance from the extremum $ i n this basis, as 5 a l; a . 
The action of the algorithm, Eq. fl53"D , on 5 a is 

5 a (r + A) = [1 - AX a ]5 a (r) , (55) 

2 It could in principle hit a saddle point and stick there, but it would have to land exactly on it and the 
space of starting points which would do so is measure zero. 
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which shrinks those perturbations in directions with positive eigenvalue and stretches per- 
turbations in directions with negative eigenvalue. This is fine if you want to find a minimum 
but it will carry you away from a saddle point. 

If on the other hand we take N steps with step length A and one step with step length 
—NA, the departure from the saddle point is 

<J a (after) = [1 + NAX a ][l - AX a ] N 8 a (before) ~ (1 + NAX a ) e - NAXa 5 a (before) , (56) 

where the last approximation holds for AA <C 1. We plot the function (1 + x)e~ x in Figure [8[ 
The algorithm reduces the deviation in every basis direction provided that the most negative 
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Figure 8: The function (1 + x)e~ x . The algorithm converges for all eigenvalues A provided 
iVAA = x satisfies -1 < (1 + x)e~ x < 1. 

eigenvalue A neg satisfies 1 + A^AA ne ge _7VAAncg > -1, which requires iVAA neg > -1.2784645. 
If we choose NA too large, the deviation in the unstable direction will alternate in sign and 
grow in magnitude; otherwise every departure from the saddle point will shrink. 

For our application we can tell if we have chosen too large a value for NA by cooling 
the configuration towards vacuum and finding the sign of TrE ■ B as the configuration slides 
out of the saddle towards the vacuum. If iVAA neg < — 1 then the sign of the excitation in 
the unstable direction switches each step, and so will the direction of approach towards the 
vacuum. If the sign switches we may be driving the system dangerously. If it does not, we 
are being too conservative and the algorithm will be inefficient. We can adjust NA as we go 
to achieve optimal stable performance. If we did not have such a clear signature for the sign 
of the unstable mode, we could still determine A neg , once we get close to the saddle point, 
by tracking the energy during cooling from near the saddle to the vacuum. After transients 
from stable modes have decayed, but before the configuration gets far from the saddle, the 
energy behaves like E sp ^ — K exp(rA neg ), and A neg can be determined by fitting. 

Also note that, if we know A neg , and also know that there is only one unstable mode, we 
can increase the efficiency with which the algorithm eliminates low frequency stable modes 
by applying extra cooling. For instance we can perform a backwards step of length 1 / A ncg 
and several small forward steps of total length 2/|A neg |, so the departure from the saddle 
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point obeys 

5 a (after) = (1 + A Q /A neg ) exp(—2A Q ,/A neg )<5 Q! (before) . (57) 

However, if we expect there to be more than one negative eigenmode then this is may make 
the other one grow unstably In general, if there are more than 1 negative eigenmodes and 
there are low lying positive eigenmodes, the algorithm becomes quite inefficient. Also note 
that the algorithm will not necessarily find a saddle point if it starts out very far from it. 
Both of these limitations are quite general limitations of saddle point finding algorithms. 
The algorithm's charm is its simplicity; all we need to be able to do is take first derivatives 
of the Hamiltonian in order to apply it, though we should try to determine A neg if we want 
to apply it well. No more coding is necessary than it takes already to search for minima. 
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